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1. Introduction and summary: naked singularities in self-similar collapse. 

In the standard picture of gravitational collapse, the implosion that results from the 
instability of the collapsing object leads to the formation of a black hole horizon prior to 
the formation of the inevitable singularity lj. However, it is known that models exist 
where the 'horizon before singularity' order is not followed, and the singularity that 
results is visible to external observers. There are different reasons why the existence 
of such naked singularities are an undesirable feature of spacetime: on a fundamental 
level, they are accompanied by Cauchy horizons leading to a breakdown in predictability 
of physical laws. On a physical level, the possibility arises that naked singularities 
may be the source of infinite (destructive) amounts of energy. The Cosmic Censorship 
Hypothesis (CCH) of Penrose seeks to guard against naked singularities. In rough terms, 
this hypothesis asserts that naked singularities cannot form in realistic gravitational 
collapse (there are rigorous mathematical formulations of the hypothesis; see for example 
|2j). Thus those models which give rise to naked singularities must be unrealistic in some 
way. Nevertheless, models admitting naked singularities provide probes of the CCH, and 
studies of naked singularities have influenced the development of exact statements of 
the hypothesis. There is also some hope that studies of spacetimes admitting naked 
singularities will shed light on how a general proof of the hypothesis might arise. Of 
course one must also keep in mind the possibility that such a model cannot be ruled out, 
and that naked singularities must be considered to be genuine astrophysical objects. 

There are at least three ways in which spacetime models admitting naked 
singularities can be ruled out as providing genuine counterexamples to the CCH. First, 
the matter model used may be considered to be inappropriate to the description of 
gravitational collapse on the smallest scale. This is understood to be the case, for 
example, in fluid models: the singularities that result are ascribed to a breakdown of 
the matter model rather than a gravitational pathology [3] . Indeed careful statements of 
the hypothesis insist that the matter model used must be such that it does not develop 
singularities in flat spacetime j2j. 

The second way is to demonstrate that the model that includes a naked singularity 
is unstable in the following way. One shows that a small perturbation of the initial data 
for the spacetime gives rise to a spacetime that does not admit a naked singularity. This 
approach has been used by Christodoulou to show that naked singularities forming in 
the self-similar collapse of a spherically symmetric massless scalar field are unstable [3]. 

Finally, one looks for instability of the Cauchy horizon associated with the naked 
singularity. In this scenario, the model admitting a naked singularity is a non-generic 
member of a class of spacetimes which instead give rise to a null singularity marking 
the end of the spacetime rather than a problematic horizon. This situation holds a the 
inner (Cauchy) horizon of charged or rotating black holes [S]. 

In this paper, we deal with a 1-parameter family of spherically symmetric self- 
similar spacetimes that admit a naked singularity: the self-similar Vaidya spacetimes. 
We will study stability of the associated Cauchy horizon. 
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The question of whether or not these spacetimes provide a serious challenge to 
the CCH can be answered immediately in the negative. The matter model is null 
dust, which always forms singularities in flat spacetime. However, there are many 
classes of spherically symmetric self-similar spacetimes admitting naked singularities 
that cannot be ruled out on this basis. By studying stability of the Cauchy horizon of 
Vaidya spacetime, we will provide a template for the study of the same issue in more 
realistic spacetimes (e.g. self-similar Lemaitre-Tolman-Bondi, perfect fluids, Einstein- 
Klein-Gordon, Einstein- SU(2)). There are also other reasons for studying perturbations 
of Vaidya spacetime. This spacetime is used to model the late stages of stellar collapse 
in which radiative emission dominates. Studying stability of the spacetime provides 
information with regard to the effectiveness of this model. 

A spacetime is said to be self-similar if it admits a homothetic Killing vector field, 
i.e. a vector field X satisfying 

£x9ab = 2g a b- 

(The choice of non-zero constant on the right hand side is arbitrary, and it should 
be noted that some authors would refer to X as defined here as a proper homothetic 
Killing vector field, or to the associated symmetry as type-1 self-similarity.) See jH] for 
an overview of the important role of self-similarity in General Relativity. 

The line element of a spherically symmetric spacetime (M, g) can always be written 
in the form 

ds 2 = 2F(v, r)e p(v ' r) dv 2 + 2e p{v ' r) dvdr + r 2 dQ 2 , (1) 

where dfl 2 is the standard line element on the unit 2-sphere. The coordinate v is an 
advanced Bondi coordinate. Taking this null coordinate to increase into the future, it 
labels past null cones of the axis r = 0. The line element above maintains its form 
under the relabelling v — > V(v). Self-similarity holds if and only if F(v,r) = G(t) and 
P(v, r) = i[)(t) for some functions G, ip and where t — -: 

ds 2 = -2G(t)e m dv 2 + 2e m dvdr + r 2 dtt 2 . (2) 

The homothetic Killing vector field isX = v-^+r-^. 

We will use the coordinates t = v/r and a rescaled radial coordinate x defined by 
r = e x . Then the line element reads 

ds 2 = e 2x {-2Ge^dt 2 + 2^(1 - 2tG)dtdx + 2te^{\ - tG)dx 2 + dQ 2 } . (3) 

A general description of spherically symmetric spacetimes modelling gravitational 
collapse has been given in [2] and jHJ- This can be done without specifying the matter 
model. Here, we restrict ourselves to the following crucial points, proven in [7j. The 
second result gives necessary and sufficient conditions for the singularity that necessarily 
forms at (v, r) = (0, 0) to be naked. 

Proposition 1 The surface t = t c constant is spacelike (respectively, null, timelike) if 
and only if t c (l — t c G{t c )) > (respectively, — 0, < 0). 
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This gives rise to the moniker 'similarity horizon' for null hypersurfaces of the form 
t = t c . 

Proposition 2 Suppose that the spacetime (M, g) with line element 

(i) satisfies the Einstein equation; 

(ii) has energy-momentum tensor satisfying the null energy condition; 

(Hi) is regular to the past of the scaling origin O = (v,r) = (0,0), where v is scaled to 
measure proper time on the regular axis {v < 0, r = 0}. 

Then there exists a future-pointing outgoing radial null geodesic with past endpoint on 
O if and only if there is a positive solution of the equation 1 — tG(t) = 0. Furthermore, 
if t\ is the smallest such positive root, then the surface t — t\ is the Cauchy horizon of 
the spacetime. ■ 



We have not defined 'regularity to the past of 0'; it suffices to note that this is a 
well-defined concept, including limiting behaviour of the metric at the regular axis and 
the absence of trapped surfaces in v < 0. Note also that we can characterise the Cauchy 
horizon as being the first similarity horizon to the future of the scaling origin. The past 
null cone of the scaling origin is a similarity horizon given by t = 0. 

In the following section, we describe the geometry of Vaidya spacetime, the subject 
of our analysis, concentrating on the self-similar case. In Section 3, we describe the 
gauge-invariant approach to perturbations of spherically symmetric spacetimes given by 
Gerlach and Sengupta PJ. We show that for odd-parity perturbations (see below), the 
matter perturbation is completely and explicitly determined by an initial data function 
(Iq. The remaining perturbation quantities are determined through a single gauge- 
invariant scalar IT, which satisifies an inhomogeneous wave equation with source term 
depending on /x . We give existence and uniqueness results, and show that, subject to 
the specification of regular initial data on a slice t — t{ G (0, t\), the function II and its 
first partial derivatives remain finite up to and at the Cauchy horizon t = t\. Thus there 
is no instability at the level of the metric or the matter perturbation. There remains the 
possibility that instability is present at the level of the conformal curvature tensor (Weyl 
tensor) - cf. the mass inflation scenario inside charged spherical black holes ^Hj • This is 
ruled out in Section 4 where we show that all the gauge and tetrad invariant perturbed 
(Newman-Penrose) Weyl scalars remain finite at the Cauchy horizon. Perturbations 
with angular mode number I = 1 require a separate (much shorter) treament, which 
is carried out in Section 5 and yields the same results. We make some concluding 
comments in Section 6. We note that the analysis of the inhomogeneous wave equation 
below follows closely our previous analysis of the minimally coupled massless scalar wave 
equation in a general spherical self-similar spacetime We use the conventions of [2|, 
and set G = c = 1. 
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2. Geometry of Vaidya Spacetime 

The Vaidya spacetime metric is the unique solution of the Einstein equations subject 
to the assumptions that (a) spacetime is spherically symmetric and (b) the energy- 
momentum tensor is that of null dust, with the dust flow vector k a normal to the 5*0(3) 
symmetry group orbits (the bar indicates a background quantity). Local conservation 
of the energy momentum tensor shows that k a is geodesic and hypersurface orthogonal, 
and so one can introduce a null coordinate v with k a = — V a v , where we assume that 
v increases into the future. The hypersurfaces v = constant are either future null cones 
of the axis {r = 0}, coresponding to expanding matter, or past null cones of the axis, 
corresponding to collapsing matter. We assume the latter (and so conform with the 
notation of the previous section). Taking as the remaining coordinates the standard 
angular coordinates (9, ip) on the unit sphere and the areal radius r, the line element 
can be shown to take the form 



where dVL 2 = dQ 2 + sin 2 9d(p 2 is the standard line element on the unit sphere. The 
energy-stress-momentum tensor is obtained from 87rt a b = pk a kb, with p = 2m'(v)/r 2 . 
Then the strong, weak and dominant energy conditions are all equivalent to m'(v) > 
for all v in the domain of m. The line element above is used to model the gravitational 
collapse of a thick shell of null dust (or a photon fluid) with the choice 



where v — > M(v) is an increasing C 1 function on [0, v\) and v\ > is arbitrary. Note 
then that spacetime is a portion of Minkowski spacetime to the past of the past null cone 
{v — 0}, and is a portion of the Schwarzschild-Kruskal space-time with mass parameter 
M(vi) to the future of the past null cone {v = v±}. The null fluid is confined to the 
region {0 < v < Vi}, and collapses from past null infinity to form a singularity at r = 0. 
The portion of the singularity in v > is future space-like, but the singular origin 
{(v,r) = (0,0)} may be visible (timelike or ingoing null), depending on the details of 
the function M at v — 0. See Figures 1 and 2f. 

In fact, we will not impose the cut off at v — v±, and will take the exterior region 
v > to be filled with null dust. That is, we take 



We wish to study the stability of the Cauchy horizon in the case when a naked singularity 
is present. In the cut-off spacetime, the portion of the Cauchy horizon that resides in 

| Note that the spacetime of Figure 2 displays a somewhat curious feature of event horizons related to 
their dependence on the global structure of spacetime: the possibility of their appearing in a region of 
spacetime whose causal past is flat. 




(4) 



r 




v < 0; 

M(v) < v < vi, 

M(vi) Vt < V, 
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r = 0, v > 




Figure 1. Conformal diagram for Vaidya-Schwarzschild collapse to a globally naked 
singularity. We use the advanced Bondi co-ordinates v and r described in the text. The 
Cauchy horizon is shown dashed, the event horizon as a double line and the apparent 
horizon as a bold curve. Note that the apparent horizon in the matter filled region 
meets the event horizon at the surface of the 'star'. Af is the past null cone of the 
scaling origin. In Figures 1-3, the matter filled region is shaded. 



r = 0, v > 




Figure 2. Conformal diagram for Vaidya-Schwarzschild collapse to a black hole. As 
in Figure 1, the apparent horizon (bold) in the matter filled region necessarily meets 
the event horizon (double line). 

the Schwarzschild-Kruskal region is (a portion of) a regular outgoing null hypersurface 
of that spacetime where no singular behaviour can be expected, unless divergence is 
mediated along the earlier portion of the horizon that resides in the matter-filled region. 
Thus we are only interested in the matter-filled region of the spacetime, and the cut-off 
is unnecessary. So our concern is the spacetime of Figure 3. 
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r = 0,v > 




Figure 3. Conformal diagram for Vaidya collapse to a globally naked singularity. The 
region to the future of v = is filled with null dust. 



This spacetime is self-similar when (and only when) M is a linear function of v. 
M(v) = Xv for A > 0. The restriction on the range of A ensures that the energy 
conditions are satisfied, and that the trivial case is avoided (A = corresponds to flat 
spacetime). Applying Proposition 2, it is then straightforward to show that the singular 
origin is naked if and only if A G (0, jg), and we will assume henceforth that A lies in 
this range. 

Introducing the coordinates t = v/r,x = logr of Section 1, the line element (0J) 
becomes 

ds 2 = e 2x [-(l-2Xt)dt 2 + 2(l-t + 2Xt 2 )dtdx + t(2-t + 2Xt 2 )dx 2 + dQ 2 ].(5) 

The Cauchy horizon is given by 

1 - v 

t = U : 



4A ' 

and the second future similarity horizon is given by 

l + v 

t = (2: =— • 

where v = Vl — 16A. There are no other future similarity horizons. The apparent 
horizon is spacelike and is located at t = t$ := 1/2A. This forms to the future of both 
future similarity horizons: t 3 > t 2 > t\. 

For the calculation of the gauge invariant perturbed Weyl curvature scalars, it will 
be essential to have an appropriate representation of the radial null directions of self- 
similar Vaidya spacetime, along with the associated null coordinates. The advanced null 
coordinate is v = tr — te x , so that v = constant describes a past null cone of the axis 
{r = 0}. The future null cones are described by u = constant, where the retarded null 
coordinate u is given by 

u = -e x '|ti -t| Al |t 2 -t\ x \ (6) 
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with 



Ai = 



v+1 



X, = 



v — 1 



1v 2v 

The region of spacetime with which we will be concerned is that bounded by past and 
future null infinity, the past null cone of the scaling origin Af = Hr and the Cauchy 
horizon H + . In the coordinates (u,v), the corresponding Lorentzian 2-space is 

M 2 = {(u,v) : -oo < u < 0,0 < v < +oo.} 

and we have the following representations (see Figure 4): 



J- = {(u,v) 
J + = {(u,v) 
H- = {{u,v) 

n + = {(u,v) 



u = — oo, < v < +oo}, 
— oo < u < 0,v = +oo}, 
— oo < u < 0, v = 0}, 
u = 0, < v < +oo}. 




Figure 4. Conformal diagram of M2- Surfaces of constant t are shown dashed, 
including the initial data surface Xj. Surfaces at infinity are shown by double lines. 

The radial null directions are ^ and J^, and we choose the following scalings: we 



will take I to be the future pointing ingoing radial null vector field given by 
r_ 2 9 

HlhL' 

where 

H = H(i) = -2A|fi - t| A2 |t 2 - t\ Xl 
and we take n to be the future pointing outgoing radial null vector field given by 



n 



d_ 

dv' 



(7) 
(8) 
(9) 



With these choices of scaling, I is an affinely parametrized null geodesic tangent vector 
field, and the Newman- Penrose normalization g a bl a n b = — 1 holds. We note also that in 
these coordinates, the line element is 

ds 2 = Hdudv + r 2 (w, v )dVt 2 . 
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3. Odd-Parity Perturbations 

3.1. The Gerlach-Sengupta formalism 

To study perturbations of this spherically symmetric spacetime, we will use the gauge 
invariant formalism introduced by Gerlach and Sengupta [H] . (We follow the presentation 
of Martin-Garcia and Gundlach [H].) This is based on the natural 2+2 splitting of a 
spherically symmetric spacetime, and a multipole decomposition that enables an efficient 
treatment of the angular dependence of the perturbation. 

The metric of a spherically symmetric space-time (Ai 4 ,g) can be written as 

ds 2 = g AB {x c )dx A dx B + r 2 (x c )^ a pdx a dx p \ (10) 

where qab is a Lorentzian metric on a 2-dimensional manifold with boundary M 2 and 
7 Q/ 3 is the standard metric on the unit 2-sphere S 2 . Capital Latin indices represent 
tensor indices on M 2 , and Greek indices are tensor indices on S 2 . r(x c ) is a scalar 
field on M 2 . 4-dimensional space-time indices will be given in lower case Latin. The 
covariant derivatives on M 4 , M 2 and S* 2 will be denoted by a semi-colon, a vertical and 
a colon respectively. €ab and e a p are covariantly constant anti-symmetric unit tensors 
with respect to g A B and 7 a/ g. We define 

VA=^ (11) 

r 

V Q = -^ + 2v A lA + 3v A v A . (12) 
Writing the energy-momentum tensor as 

t ab dx a dx b = t AB (x c )dx A dx B + Q(x c )r 2 ^ a(3 dx a dx p } (13) 
the Einstein equations of the spherically symmetric background read 

Gab = - 2(v A \b + v A v B ) + V g A B = Sirt A B (14) 

X -Gl = - K + v A lA + V \ A = : 8vrQ, (15) 

where G^ = 7 Q ^G ai g and TZ is the Gaussian curvature of M 2 . 

Spherical symmetry of the background allows us to expand the perturbed metric 
tensor in terms of spherical harmonics. Writing Y = Y™ and suppressing the indices 
l,m throughout, we have the following bases for scalar, vector and tensor harmonics 
respectively: {Y}, {Y a := Y. a , S a := e^Yg} and {Y>y a p, Z a/3 := Y a:/3 + l -^-Y^ a/3 , S a:/3 + 
Sp-.a}- These are further classified depending on the transformation properties under 
spatial inversion x —>■ —x on the unit sphere: a spherical harmonic with index I is called 
even if it transforms as (—1)' and is called odd if it transforms as (— 1)' +1 . In the bases 
above, Y,Y a and Z a p are even and S a ,S( a: fl are odd. We note that in the analysis 
below, the multipole index I only appears in the combination (/ — + 2) and so we 
define t= (/-l)(/ + 2). 

The perturbation 5g a b of the metric tensor can then be decomposed as 

SgAB = h AB Y, (16) 
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SgAfi =h E A Y :(3 + h° A S Pl (17) 
8g aP = r 2 K laP Y + r 2 GZ a/3 + 2hS {a:f3) . (18) 

The superscripts E,0 stand for even and odd respectively. Note that h AB , {h A ,h A } 
and {K, G, h} are respectively a 2-tensor, vectors and scalars on M 2 . A similar 
decomposition of the perturbation of the stress-energy tensor is made: 

5t AB = At AB Y, (19) 
5t Af3 = AtfY.p + At°Sp, (20) 
6t a/3 = r 2 At 3 lal3 Y + r 2 At 2 Z a(3 + 2AtS {a:/3 y (21) 

In this case, At^B, {At A , At A } and {At 3 , At 2 , At} are respectively a 2-tensor, vectors 
and scalars on M 2 . 

A complete set of gauge invariant variables is produced as follows. An infmitesmal 
co-ordinate transformation on the background is generated by a vector field £. Again, we 
can decompose into even and odd harmonics and consider separately the transformations 
generated by the 1-form fields 

£ E = £ A ( x c )Ydx A + i E {x c )Y.. a dx a , (22) 
f = i°S a dx a . (23) 

From the transformed versions of the metric perturbations, one can construct 
combinations which are independent of the coefficients of £. These combinations are 
then gauge invariant. As we will only be concerned with the odd parity sector, we give 
only these terms. The entire odd parity metric perturbation is captured by the gauge 
invariant co- vector field 



k A = h A — h\ A + 2hvA, 

and the gauge invariant matter perturbation is described by 

L A = At° A - Qh c , 
L = At- Qh. 

The linearized Einstein equations read 



1< A- 



k A 



— r 



kc 



- ~^k A 



(24) 

(25) 
(26) 

(27) 



k A 



\ A = 16vrL. (28) 

The latter equation follows from the former and the linearized conservation equation 

(r 2 L A y A = £L. (29) 

As noted by Gerlach and Sengupta 9 3 , the vector equation (|2Tj) is equivalent to the 
single scalar equation 

1 



(r 4 n)' 



-m = -167T6 



AB t 

L A\B, 



(30) 
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where 

U:=e AB (r- 2 k A ) lB (31) 

is a gauge invariant scalar. k A is recovered from II and L A using 

ik A = 16vrr 2 L A - e AB (r A U) lB . (32) 

The quantity II is not only a gauge invariant scalar that, along with L A , completely 
determines the metric perturbation, but as shown in [T2], has the tetrad and gauge 
invariant interpretation of being the perturbation of the Coulomb component ^2 of the 
background Weyl tensor. 

In order to close the system of perturbation equations, an equation of state must 
be given for the perturbed spacetime. With this addition, the equations (J2HJ an d PHjl 
completely determine the perturbation. 

An important point to note is that the formalism described above is incomplete for 
1 = 1. (There is of course no odd parity I = perturbation.) For I = 1, h is not defined, 
being a coefficient of zero, and so should be considered to be zero. Thus the gauge 
invariants cannot be constructed. However it is convenient to use the same variables 
()24 )) -([26 )) for all values of I. For I = 1, these variables are only partially gauge invariant, 
and so gauge- fixing is required. We defer treatment of the case I = 1(£ = 0) to Section 
5, and so assume until then that I > 2. 

3.2. The matter perturbation 

We assume that the non-vacuum portion of the perturbed spacetime is filled with null 
dust. This allows us to write 

tab = {p + Sp)(k a + 5k a )(k b + Sk b ), 

where the barred terms refer to background quantities, so that k a = — V a v. Retaining 
only first order terms, we define 

St ab = Spk a k b + p(k a 5k b + k b 5k a ). 

Comparing with ([19| ) -([2i p and recalling that we are considering only odd parity 
perturbations, we can determine the gauge invariant matter perturbation using (|25j) 
and (J2HJ): we find L = and 

L A = p(x A )pk A) (33) 

for some (first order) scalar p. The evolution of p is controlled by (j29|) . In the self-similar 
coordinates (t, x) we find 

tp,t ~ P,x = 0, (34) 

which yields 

p = p (te x ) = p {v) 

for some arbitrary different iable function pq. Thus the matter perturbation is completely 
determined by the specification of the function p Q on an initial slice of the form 

t = u e (o,ti). 
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3.3. The master equation 

Having completely specified the matter perturbation in terms of an initial data function, 
we now turn to the only remaining perturbation equation ()3()j) . We will refer to this 
as the master equation and will work in the coordinates (t, x) of Section 2. It is worth 
repeating that we deal only with the region < t < t\ on which t is a time coordinate. 
The analysis below follows closely that of jH], and where possible, we will quote results 
from this paper rather than repeating very similar proofs. We reiterate that I >2. 

For fixed k G R, we define 0(t, x) = e KX U = r K U. Then, using (|3U|) and the line 
element (J5J), we find that satisfies the inhomogeneous wave equation 

a<f>, a + + 70,^ + («' + (6 - 2«)/3)£ t + (/?' + (6 - 2k) 7 )0^ 



+ ((4 - k)$ + (8 - 6k + k 2 ) 7 - te*)4> = 8\e {K - 3)x fi {te x ), (35) 
where if) = and 

a = -2t(l-tG), (36) 

P = 1 - 2tG, (37) 

7 = 2G, (38) 

G = i - At. (39) 



Modulo the specification of the functional form of G, the left hand side of (|3~5j). along 
with (J36|) -(|38| ) . gives the general form of the left hand side of the master equation for 
the line element (j2J). 

We fix an initial data surface for (j33|) given by £; = {(i, x) : t — t^x G M} with 
ti G (0, ti). Our principal concern is how and various of its derivatives representing 
gauge invariant curvature scalars behave in the approach to the Cauchy horizon, subject 
to initial regularity conditions imposed at Sj. Noting that we may write 

a{t) = -2Xt(tx -t)(t 2 -t) 

and so a(ti) = 0, we see that this question is rendered nontrivial by virtue of the fact 
that the Cauchy horizon is a singular hypersurface for the equation ()35|): the spacelike 
surfaces t =constant become characteristic (null) in the limit t — > t\. Prior to the 
Cauchy horizon, the evolution of proceeds smoothly: 

Theorem 1 Let f,g,p,o G G^°(R) and let k el Then there exists a unique solution 
<f) G C°°([ti,ti) x R, R) of the initial value problem consisting of the equation A 3b}) and 
the initial data 

0k = /, 0,tk = g- 

Furthermore, the solution satisfies <p\t=t G Gq°(R) for all to G [£,,£].). ■ 
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The proof of this result is standard and is most easily obtained by rewriting the 
equation (J33j) as a first order symmetric hyperbolic system for 



<P 



a<f) 7t + (3(j). 



\ 



',x 



\ 



J 



(40) 



See for example Chapter 12 of [T^]. It is convenient to rescale the time coordinate by 
defining 

rl ds 



a[s) 



(41) 



Then r is an analytic function of t on ti), r(tj) = and limi_> tl r = +oo. The master 
equation (J33j) can be written in first order symmetric hyperbolic form 

tT + A0 iX + Bip = j 

where A, B are smooth, bounded matrix functions of r on [0, oo) and A is symmetric 
with real distinct eigenvalues. The source term j is given by 



V 



o \ 

]\a(t)e^-^ x ^ (te x ) 

I 



We wish to analyse the behaviour of the solutions described by Theorem 1, and so 
we assume until indicated otherwise that the hypotheses of this theorem hold. Then 
JG C °°(Mi) x B 

We define 



EAt) 



dx. 



The growth of this energy norm is described by the following corollary, again a 
standard result. We use the notation ||a|| 2 = (a, a) (Euclidean inner product) and 
||a||| = j R \\a\\ 2 dx (L 2 norm). The terms Cq, C\, ... represent possibly different constants 
that depend only on the metric function G and the angular mode number I. 

Corollary 1 £' 1 [0](r) is differentiate on [0, oo) and satisfies 



E 1 [<p}(r)<e B ^(E 1 [<f ) }(0) + 



WWW, 



(42) 



where B Q = sup r>0 \I — 2B\ < +oo. Consequently, 



\<p(t,x)\ 2 dx < e BoT (Ei[0](O) 



,(t,x)\ 2 dx < e BoT (E 



;t(t,x)\ 2 dx KC^iE^m 



\W)\W, 



{a)\\lcb), 
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The bounds on the L 2 norm of (J) and its derivatives come straight from the definition 
of Ei(t): the third requires the use of Minkowski's inequality and incorporates the 
exponential growth of a -1 as r — > oo. As in jH], the growth of these norms in the 
approach to the Cauchy horizon r — > oo is analysed using a second energy integral. 

Let 

D{t) =£+( K - A)(3'(t) - (k 2 - 6k + 8) 7 (t), 
and for an arbitrary positive, real-valued differentiable function K(t) define 

E 2 [0, fi ] (t) = [ -a<f) 2 t + j<t>l + D<P 2 + Ke 2{K - 3)x fi 2 dx. (43) 
Jr 

Lemma 1 Let k G [0,4]. Then for all t G D(t) > and D'(t) < 0. 

Proof: From the definitions (|37p and (|38|1 we obtain 

D(t) = £ - (k 2 - 5k + 4) + 2\t(K 2 - 4k). 

Thus D'(t) = 2Ak(k - 4) < 0. Then for t < t x , 

D(t) > D{tx) = £ - 4 + (5 - 8A*i)« - (1 - 2Ati)K 2 
> (5-8At!)K-(l-2Ati)K 2 , 

where the second inequality uses / > 2. This last expression, considered as a quadratic 
function of k, is non-negative for k G [0, «*] where 

5 - BXti 2 

k^ = — = 4 H > 4. 

1 - 2Ati l + i/ 

This yields D(t) > on the range indicated. ■ 



Lemma 2 Let k G [0, Ki) where 



1. ^ . 2 , 1 + 32A - 32A 2 - v 
ki = -(5 - 4A + At 2 ) - 



2 V Ll 16A 

T/ien t/iere exzsts t c = t c (n) G [tj,ti) ; a constant C and a choice of K such that 

IE2 
dt 



E 2 [(j)^o}(t) > and ^ < C ^2 /or all t G [t c ,t 



Proof: Noting that k x < 4, we see that Lemma 1 applies and so non-negativity of E 2 is 
immediate. Write fi = e^ K ~ 3 ^ x fXo- E 2 (t) is a smooth function of t, and smoothness of the 
solution cf) and of Hq allow differentiation under the integral sign. The resulting integral 
is simplified in three steps: (i) integration by parts of the term §,x<\>,xt and the removal 
of a boundary term - permitted as <fi has compact support on each slice t = constant; 
(ii) removal of the term with cf) M by application of the equation ()35|) ; (iii) removal of a 
total derivative containing <pj4>,xt- This results in 
dE 2 



dt 



[ [(a' + 2(6 - 2k))m 2 + 2{f3' + (6 - 2k) 7 )w t w , + j'u% + D'u 2 
Jr 

— I6\jlu tt + 2Kfifi t + K'ji 2 ] dx. 
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In the next round of simplifications, we apply Lemma 1 (D' < 0), the Cauchy-Schwarz 
inequality in the form 

/ 2jlu t t dx < / fi 2 + u\ dx 
Jr Jr 

and the equation satisfied by fi (obtained from (JMj)) 
tjlj ~ fi,x + — 3)/i = 0. 

These yield 

dE f 

'' ' 2 < / [(«' + 2(6 - 2k) + 8X)u 2 t + 2{f3' + (6 - 2 K )j)u, t u, x + -f'u 2 x 



dt 



+ (^'_ 2 ( K -3)^ + 8A)/i 2 

v 



I dx. 



dx 



For any constant C > 0, define Ir by / = CIe 2 + Ir, where Ie 2 is defined so that 
E2 = L Ie 2 dx. That is, 

I R = (a' + 2(6 - 2k) + 8A)«J + 2(/3' + (6 - 2k)7)u t u >a; + YuJ. - CDm 2 
+ (Hf - (C + ~(k - 3))if + 8A)/t 2 . 

The Lemma is proven by showing that there is a choice Co > of C and a value 
t c E [ti,tx) for which I R < on [i c , tj.). 

For any choice of C, there is a positive differentiable function K defined on 
for which the coefficient of ft 2 in is negative for all k in the range specified. Making 
this choice and applying Lemma 1, we obtain 



Ir < («' +Ca + 2(6 - 2k) + 8\)u 2 t + 2{f3' + (6 - 2K)j)u t u >x + (7 - C^u 2 , 
= : a(t)u 2 t + b(t)u )t u :X + c(t)u 2 x . 

Consider next the quadratic form 

Q(X, Y; t) = a(t)X 2 + b(t)XY + c(t)Y 2 . 

We find 

a{t x ) = 2(4/t-5 + 4A + At 2 ), 
b{tx) = -i(2 K -5-At 2 ), 

c{tx) = -2(A+£). 

We note that c(ii) < for any C > and A G (0, 1/16). The term a(ti) is negative due to 
the assumed bound on k. Then the discriminant is also negative if C is chosen sufficiently 
large. So Q(X,Y;ti) is negative definite. Then by continuity of the coefficients a,b,c, 
the quadratic form Q(X, Y; t) = a(t)X 2 + b(t)XY + c(t)Y 2 is negative definite for all t 
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sufficiently close to t\. That is, there exists some t c < t\ such that Q(X,Y;t) < for 
all and t G [t c , ti) with equality holding iff X — Y — 0. We note however that 

the value of t c will depend on k, with t c — > t\ This however does not affect 

the proof, which is now completed. ■ 



Remark 2.1 By minimising K\{\) for A G (0,1/16), we could restate Lemma 2 with 
the simpler requirement k G [0, |). 

We can now give our first main result. 

Theorem 2 Let <fi be a solution of h3ty that is subject to the hypotheses of Theorem 1 
and Lemma 2. Then the energy E 2 [(f), £*o](*) of the solution satisfies the a priori bound 

e 2 [<I>, fi ](t) < Ci^ko) + c 2 jm, t g Mi) 

where 

JM= [ e 2 ^ x ^(Ue x )dx. 
Jr 

Proof: We point out first how to convert the bounds of Corollary 1 to a priori bounds. 
To do this, we exploit the self-similar nature of the solution /i of the matter perturbation 
equation (JMj) . We have 

\\j]\l(r)=64a\t) [ e 2 ^ x ^{te x ) dx. 

A change of variable in the integral yields 



a(ti) ) \t 

Then making the appropriate change of variables using ([41)1 . we obtain 

=: h(t)J K [iio\ < C J K [fi }. 

From Corollary 1, we can write 

Eh[<j>,iM>)(t) < d(t)(£7i[0](O) + C 2 JM), (44) 

where d(t) is a smooth positive function of t that diverges in the limit t — > t\. However, 
d(t c ) is finite, where t c is the value of t identified in Lemma 2. From time t c onwards, 
E 2 obeys the differential inequality of this lemma, which may be integrated to yield 

E 2 (t) <e Co ^E 2 (t c ), te[t c ,h). (45) 

Combining (jUJ) and yields the desired result. ■ 
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Theorem 3 Let <fi be a solution of \35\) that is subject to the hypotheses of Theorem 1 
and Lemma 2. Then cf) is uniformly bounded on [ti, £i) x R: there exist positive constants 
C\ and C*2 such that 

\<f>(t,x)\ < CtE^iO) + C 2 J K [pL ], teMO.zeR. (46) 

Proof: Theorem 2 provides an a priori bound for the if 1,2 norm of (f>: for all t G [ti, ti), 

[ 2 + </> 2 dx < C X E X [4>]{$) + C 2 JM- 
Jm 

The pointwise bound arises immediately on application of the Sobolev inequality for 

v G C£°(R): 

\v(x)\ 2 <±yjv(y)\ 2 + \v'(y)\ 2 dy\. 
See p. 1057 of [TJ] for a proof of this inequality. ■ 

Remark 3.1 Note that by differentiating with respect to x, we can obtain results 
similar to Theorems 1-3 for any spatial derivative of <fi. The bounding term in the 
inequalities corresponding to (|4~E|) will involve sums of terms of the form £'i[|^](0) and 

Remark 3.2 Of principal concern in this paper is the behaviour of the field II and 
those of its derivatives representing the perturbed Weyl curvature scalars. Theorem 3 
shows that II is bounded in the limit as the Cauchy horizon is approached (t — > ti). 
However this does not imply that lim t _* tl Tl(t, x) exists for any x G R. We will show 
now that this is in fact the case; indeed we can show that II| t=tl G C°°(R). In [Hj, 
the corresponding limit function was erroneously assumed to exist on the basis of the 
equivalent to Theorem 3. This assumption can be shown to be true by applying the 
argument below to that paper, and so does not affect the results of that paper. 

In order to get from the bound of Theorem 3 to the existence of the limit 
lim t ^ tl (f)(t, x), we need some control over the time derivative of as the Cauchy horizon 
is approached. This is provided by the following lemma, which relies on treating (f33|) 
as a first order transport equation for <p )t . 

Lemma 3 Let </> be a solution of hSfy) that is subject to the hypotheses of Theorem 1 and 
Lemma 2. Then <p t is uniformly bounded on [tj,ti) x R: there exist positive constants 
Co, Ci, ... such that 

|0, t (t, x)\ < CoE^Mp) + C 1 E 1 [<f )tX ](0) + C^^KO) 

+ c 3 JM + C 4 JM + C s J K \p®> t e [U, h),xe R. 



Proof: Define \ = <p t t- Then (|35|) can be written 

a X ,t + Wx,x + («' + (6 - 2k)(3) X = fit, x), 



(47) 
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where the right hand side depends linearly on fi and the zeroth, first and second spatial 
derivatives of 0. Then / is smooth and has compact support on each slice t =constant. 
If we write ()33|) as 

L[<j>] = ao(t,x)fx (te x ), 

then differentiation with respect to x shows that )X and (j) yXX satisfy equations with 
identical first and second order derivative coefficients: 

L[4>,x\ = b (t, x)fio(te x ) + bx(t, x)fi' (te x ), 

L[(j),xx\ = c (t,x)fi (te x ) + ci(t,x)n' (te x ) + c 2 {t ) x)^(te x ). 

We can therefore apply Theorem 3 to <p yX and <p )XX : the only difference in the result will 
be that the bounding terms will depend also on the L 2 norms of the first and second 
derivatives of <fi and Hq at t = ti. Then by linearity, we can bound / by a similar a 
priori term. The bound for u then arises by straightforward integration of the first order 
transport equation (JTTj). See Theorem 6 of [H]. ■ 



Theorem 4 Let <fi be a solution of h3ty that is subject to the hypotheses of Theorem 1 
and Lemma 2. Then (fi n + : = lim^ tl <p(t, •) G C°°(R) and satisfies the bound 

\<Ph+(x)\ < Ci^i[0](O) + C 2 JM, xeR. 

Proof: Fix x G M. and consider a sequence {t ( - n ^}^L 1 C [U,ti) that converges to t±. For 
all m, n > 1, we can apply the mean value theorem to get 

|0(t (m) ,x) -0(t (n) ,x)| = \Mt*,x)\\t (m) ~t {n) \ (48) 

for some between t^ and t^ n \ Using the bound of Lemma 3, we see that (j)(t^ m \x) is 
a Cauchy sequence of real numbers, and so for each x G 1R, lim^^ <p(t, x) exists. Hence 
4> n + is defined. We can apply an analogous argument to all the spatial derivatives of (j>. 
It remains to show that 

^-{<Ph+} = Ihn (j), x (t,x). 
ax *-»ti 

(Again, an analogous argument will apply to higher spatial derivatives.) But this follows 
by uniform convergence of the sequence of functions (p(t^ n \x) to (fin+i which in turn 
follows from ()48|) and the uniform bound of Lemma 3. To obtain the bound in the 
statement of the theorem, we take the limit t — > £i of the corresponding bound in 
Theorem 3. This is permitted as the bounding term is independent of t. ■ 

We conclude this section by extending the results of Theorems 1-4 to the case where 
the initial data lie in appropriate Sobolev spaces. This is important as it will allow the 
perturbation to be non-zero at the axis r = 0. The results so far relate to the case where 
the initial data and the corresponding solutions are supported away from r = 0. This 
is an undesirable feature, as we would ideally like to consider a perturbation that arises 
from data imposed on a globally regular initial data slice of the space-time: our t = t { 
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slice is singular at the origin. Such a regular slice would intersect the past null cone, 
and we should certainly consider the case where the support of the initial data also does 
so. We can deal with such data by taking the limit of a sequence of test function (Cq°) 
data in an appropriate Sobolev space. 

Theorem 5 Let k G [0, Ki) and let <p — e Ka n and define jl(x) = e*-" -3 ^ n$(x) = 

(i) Let f G H 1,2 (M), g G L 2 (R), fi G L 2 (R). Then there exists a unique solution 
4> G C([ti, ti), iJ 1,2 (R)) of the initial value problem consisting of Iff 5)) with initial 
data 0|si = /, 0,t|sj = 9- The solution satisfies the a priori bound 

\<f>{t,x)\ < C/ o £i[0](O) + C 2 J K [fl}, t G [ti, ti),x G R. 

(ii) Let f G H 3,2 (M,), g G H 2 ' 2 (R), fi G H 2 ' 2 (R). Then there exists a unique solution 
4> G C([£j, ti], H 1,2 (M)) of the initial value problem consisting of < 1,?,5|) with initial 
data = /, 0,t|si = 9- The solution satisfies the a priori bound 

\<f>(t,x)\ < Co^i[0](O) +C 2 J K [fi], t G Mi], are R, 

and its time derivative satisfies 

\Mt>*)\ < Co^i[0](O) +C 1 E 1 [<f> >x ] (0) +C 2 E 1 [<l>^](0) 

+ C 3 J K \fi] + C A J K [jl'} + C 5 J K \jl"}, t e [U, ti],x G R. 

Proof: We give just the outline of the proof, which is nearly identical to Theorems 5 
and 7 of [8 a and which relies on a standard PDE technique. For part (i), we consider 
sequences of test functions {/( n )}^L , {g( n )}^L , {/i( n )}^ and apply Theorems 1-3 to 
obtain a sequence of smooth solutions {0( n )}^o c C°°([ti, ti) x R) satisfying the bounds 
of Theorems 2 and 3 above. By exploiting linearity of the equation and the fact that 
Co°(R) is dense in the Banach spaces L 2 and H k ' 2 for k = 1, 2, 3 . . ., we can legitimately 
take the limit of relevant inequalities to prove the stated results. The proof of part 
(ii) is similar, but higher order Sobolev spaces must be invoked due to the form of the 
inequality of Lemma 2. ■ 

4. Gauge Invariant Curvature Scalars 

As seen above, the odd parity perturbation of Vaidya spacetime is completely described 
by the vector La and the scalar n. The latter quantity plays a dual role: on the one 
hand it is a potential for the gauge invariant metric perturbation kA (see (I32j) ) and on 
the other, is the gauge invariant perturbation of the Coulomb component ^2 of the Weyl 
tensor [T2]. For at least two reasons, it is desirable to have a full set of gauge invariant 
scalars that describes the perturbation of the Weyl tensor. 

First, one prefers scalars as these avoid the problems presented by an inappropriate 
choice of coordinates. The components of a non-zero rank tensor may blow-up, with 
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the blow-up inadvertently ascribed to singular behaviour rather than the wrong choice 
of coordinates. 

Second, the metric and matter perturbations alone do not capture the whole 
physical picture (neither of course does the Weyl tensor alone). Perhaps the best 
example of this is the case of perturbations impinging on the inner (Cauchy) horizon 
of a charged or spinning black hole. Here the metric perturbation remains continuous 
(but not differentiable) and the Weyl curvature blows up, a scenario described as mass 
inflation This has been described rigourously in [TH] . 

As shown in [T2], the perturbed Weyl scalars can be defined in a gauge invariant 
manner in the case of odd parity perturbations. Using a null tetrad {I, n, m, fh*} where 
the asterisk represents complex conjugation and following the notation of ^H] for the 
Weyl scalars *q, . . . , *4, the perturbed Weyl scalars are given by 

Qo 



5* 
5*i 

5* 2 
5* 3 



— l A l B k AlB , 

Qi 



r 

Ql 

r 

5*.= ^1 



(r 2 n), A / 



-,kd A 



(r 2 U) ]A n A - -k A n J 



(49) 

(50) 
(51) 
(52) 



2r 2 



A^B h 



n n 



(53) 

where the functions Qi,i = 0,1,2 depend only on the angular coordinates. In this 
definition, we restrict to the preferred sets of null tetrads defined on the spherically 
symmetric background for which I and n are the principal null directions of the 
background. Then there remains a scaling freedom in these definitions. Under the 
spin-boost transformation 

(I, ft) -> (Al, A" 1 ft) 

we find 

5*; -> A 3 " 4 6% 

for i — 0, . . . 4. Thus we cannot ascribe direct physical significance to the values of 5*j 
(except for i = 2). However the terms 

r /2 , (54) 

(55) 

1 1/2 (56) 

are fully invariant perturbation scalars: i.e. they are first order scalars which are both 
gauge and tetrad invariant. 



SP-i 


= |5* 5* 




= 5* 2 , 


SPi 


= |5* x 5* 



4-1. The master equation in null coordinates 

It is straightforward to calculate 5* ,...5* 4 in the coordinates (t,x) of Section 3. 
However it is less straightforward to determine whether or not the resulting quantities 
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are finite at the Cauchy horizon t = t\. we encounter terms involving products of 
negative powers of the term H(t) of Section 2 with first and second time (t) derivatives of 
IT. While we know that the first time derivative of II is bounded, we have no information 
regarding the behaviour of the second time derivative at the Cauchy horizon. It turns 
out that we can circumvent this problem by calculating the scalars (j49j) - (j53)l in the null 
coordinates (u, v) of Section 2. This approach also necessitates rewriting the equation 
in these null coordinates. The result of this is as follows. Let 

* = r 3 n = e 3x n. 

Then we find 

tf,™ + V(u,v)V = F, (57) 



where 



V= - —(a(3' + 1 -(3 2 -to), 
Auv 

01 

F= -2\—no(v). 

uv 

In these coordinates, we find that 

4 2 — a' 

<J *o = -T7j(*«« + -^— (58) 

for which the problem mentioned above remains. However the following relabelling of 
the null cones resolves the problem. 
Define 

X = \u\^ = \ti -t\\t 2 -t|^exp(^), Y = v = te x . (59) 

We note that X = at the Cauchy horizon t = t\. Then ()57j) - that is, the master 
equation - takes the form 



where 



and 



y iXY + V(X,Y)V = F(X,Y), (60) 



V(X,Y) = -^A 1 A|t 2 -t|^exp(-(i^)x) 



1 I Q 

F(X, Y) = 4AiA 2 |t 2 - t\& exp(-(^-)x)^ (y). 

Note that for any fixed Y = v > 0, V and F are analytic functions of X at X = 0. 
Moreover, there exists C > and sequences of smooth functions {Ki}^o> {^n}^=o sucn 
that 

oo oo 

V(X,Y) = J2V n (Y)X n , F{X,Y) = Y,F n {Y)X n , \X\ < C. (61) 

n=0 n=0 

Combining the high degree of regularity of the coefficients of (jfiOj) with the results of 
Section 3 yields the following. 
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Theorem 6 Let f } g,^ G C^°(R) and let II G <7 00 ([t i ,ti) xl)fl = t ± } x R) 5e 

the unique solution of the initial value problem consisting of \35)) with initial data 

n| Si = /, n t | Ei = g. 

Let X, Y be as defined in Then there exists X > such that \1> = e 3x n satisfies 

* G w/iere 

ft = {(X, F) : < X < X , Y > 0}. 

Thus all X— and Y— derivatives of ^ are finite at the Cauchy horizon X = 0. 

Proof: Consider the characteristic rectangle 

R = {(X, F) : < X < X , Y < Y < Y 1 ] 

where Xo > and < Yq < Yi. Applying Theorem 1 and Theorem 4 with k = and 
noting that the coordinate transformation (t, x) — > (u, v) is a homeomorphism on R for 
sufficiently small X , we see that ^ G C°(R, M). Then rewriting ()60|) as 

* lX y = -y* + F=:Q(X,r), 

we have Q G C°(R, M). For any F G (Y , FJ, we can then integrate to obtain 

We can choose Yq to be small enough so that the ingoing null ray Y = Yq lies outside 
the support of \& (see Figure 5). Then f]|(X, Y ) = 0, and 

— (X,Y) = Q(X,Z)dZ, 

giving V t x G C°(i?, R). A similar argument yields ^f iY G C°(i?,R), and so * G 
C 1 (i?, R). Continuing this argument inductively, we obtain \I/ G C°°(R, R) for all Fo 
taken sufficiently small and all Y\ > 0. ■ 

Remark 6.1 It is reasonable to ask why we have not simply stated the entire problem in 
the coordinates (X, F) and deduced finiteness of ^ - and hence II - at the Cauchy horizon 
by writing down a very simple existence and uniqueness theorem for the characteristic 
initial value problem consisting of (|5Uj) with characteristic data ^\x=x ^\y=y G 
C°°(0, +oo). The answer to this is that this formulation of the problem assumes that the 
field II is a smooth function of the retarded null coordinate X at the Cauchy horizon. 
But the question of whether or not this is a valid assumption is the very point that we 
are attempting to address here, with respect to finite initial Cauchy data posed on a 
hypersurface that precedes the Cauchy horizon. 
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Figure 5. Spread of the support of the field (shaded), with initial data with compact 
support So C E». There are ingoing null rays from Ej to H + which lie outside the 
support of the field. 



4-2. The perturbed Weyl scalar s 

There is a further advantage of the coordinates (X,Y). Not only does the master 
equation assume the simple form ([HO]), but we find that the perturbed Weyl scalars 
assume a very simple form when expressed in these coordinates. Using f!58|) and (JoTJj) . 
we find 



M = -^X^(*,xx + (1 - Y a ') X ^,x)e- x Q . (62) 



Using flSJl, flMl), and ((HEJ), we have 

if- 2 X— = —\t 2 - t\-~ exp(2(— —)x), 

(1 - |a')X^ = AxA(2t 2 - i x - 3*)(* 2 - exp(-(^)x). 

Thus 5^/o is finite at the Cauchy horizon X = 0(t = This also holds for the other 
perturbed Weyl curvature scalars. Finiteness is immediate for 5^2, which is essentially 
II. The other gravitational radiation scalar, representing outgoing radiation, is found 
to be 

4 = {(^ iyy +(l-4At)e-^ i y-4Ae- 2 ^)e :E -4A(/i[ ) (r)-2Ate- :c /i (r))}e- 2a; Q*.(63) 

It is immediate from Theorem 6 that this is finite at the Cauchy horizon. For the two 
remaining scalars, we find 



with 



5^f 1 = ||_(£-4)iJ- 1 X^ 1 e^,x + (^ + 4)^|e- 3x ^ 



H -i X ^ = -—\t 2 - t\~^ exp((£-p]» 



(64) 
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while 

5^3 = j(* + 4)* y - - 4)*(1 - 2Xt)Y~ 1 ^ + 647rpe x, /i (V)} e" 2a; ^i.(65) 

Thus these scalars are also finite at the Cauchy horizon. Thus we have proven the 
following. 

Theorem 7 Subject to the hypotheses of Theorem 6, the perturbed Weyl scalars \5(A) - 
\5 l J\) . calculated with respect to the null vectors {l,n} defined in and (0J) ; are finite 
at all points of {(t,x) : ti < t < t\,x G K}. /n particular, the perturbed Weyl scalars 
are finite at the Cauchy horizon 7i + . ■ 

Remark 7.1 Having found a null tetrad in which all the perturbed Weyl scalars are 
finite at the Cauchy horizon, it is clear that the fully invariant scalars 5 Pi, i = 0, ±1 are 
also finite thereat. A difficulty of interpretation of the 5^5/ i would only arise if one found 
(say) vanishing of 5^/ and divergence of at the Cauchy horizon. In such a case, 
recourse to the calculation of <5P_i would be essential. 

As in the previous section, we wish to extend the present results to initial data 
of the type considered in Theorem 5. As in the proof of that theorem, the principal 
requirement is the existence of a priori bounds for the relevant quantities. 

Theorem 8 Subject to the hypotheses of Theorem 6, the following a priori bounds hold: 
for each i 6 {0, . . . , A], there exist constants Cq, C[ l \ . . . , Cf 1 that depend only on A 
and k such that for all (t, x) G [ti, t\] x M ; 

|e K3 W 4 (t,x)| < Cf £i[0](O) + C?E 1 [ ( f> >x }(0) + C?E 1 [ ( f> >xx }(0) 

+ Cf J K [H + Cf J, [// ] + Cf } J K . (66) 

Proof: We consider first the case k = 0. For % = 2, the bound ()66)1 is immediate from 
the definition (j51j) and from Theorems 2 and 3 (where we let k = in those theorems). 
A straightforward calculation shows that 

«j* 1 = [a(t)n, t + 6(t)n !e + c(t)ii]^, 

where here and in the rest of this proof, a,b,c, . . . represent functions of t that are smooth 
and uniformly bounded on [ti, ti], and which may change from line to line. Similarly, 
we find 

= [o(*)n j( + b(t)U x + c(t)U + I287r\e- 3x n (te x )] ^. 

The result follows by application of Theorems 2 and 3 to II and analogous results for 
U jX (see Remark 3.1) and by application of Lemma 3 which provides the bound for 11^. 
Again, we take k = in these results. 
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To obtain bounds on and 6^4, we exploit the form (|6U|) that the master equation 
takes in the null coordinates (X, Y). Integrating, and using the fact that \& is identically 
zero for sufficiently small values of Y, we can write 

V iX {X,Y)= [ -V(X,Z)V(X,Z) + F(X,Z)dZ. 
Jo 

Differentiating under the integral sign (which is permitted by smoothness) gives 

*,xi= / -V X V + V^ tX + F >x dZ. 
Jo 

This can be written as 

f Y 1 — v 

*n = / [a(t)U t + b(t)U x + c(t)II + d(t)e~ 3x ^{Z)\ exp(2(— -)x) dZ, 
Jo 1 + v 

where in the integral it is understood that t = t(X, Z) and x = x(X, Z). The term in 

the integrand in square brackets can be bounded by an a priori term of the form in the 

statement of the theorem; we use M to represent such a term. Then 

Wxxl < M I Z 2 ^) dZ, 







wher we have used the definition Y = te x and absorbed a function of type a(t) into M. 
Evaluating the integral, we see that 

\H- 2 X^e- x m :X x\ < a(t)M = M. 
In a similar manner, we can show that 

|(1 _ ^ a ')H- 2 X-^e~ x ^ iX \ < M. 

Hence by f!62|) . the theorem is proven for i = 0. The case i = 4 is similar (but more 
straightforward) . 

For values of k with k ^ 0, we simply point out that by virtue of the self-similar and 
linear nature of the equation (|H5]l and the linearity of the in II and fi, an identical 
argument to that above applies. ■ 

We conclude by writing down the result describing bounds on the perturbed Weyl 
scalars that ensues by considering data of the form dealt with in Theorem 5. We omit 
the proof, noting that this proceeds in exactly the way described in the summary proof 
of Theorem 5: we apply Theorem 8 to a sequence of solutions generated by data in 
Cq°. Then the hypothesis that the limit of the data exists and lies in (some) H n ' 2 and 
the existence of the a priori bounds in Theorem 8 ensures the existence of the limits of 
those bounds. 

Theorem 9 Let k e [0, Ki) and let <f> = e KX U and define jl(x) = e*-" -3 ^ Ho(x) = 
e^ x n(tie x ). Let f E H 3 ' 2 (R), g E H 2 ' 2 (R), ft E H 2 ' 2 (R). Then the perturbed Weyl 
scalars f^ )- < f33|) calculated with respect to the unique solution (f> E C([ti,tx], /f 1,2 (R)) of 
the initial value problem consisting of KSb^l with initial data = / ; 0,t|sj = 9 satisfy 
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the a priori bounds ib'b]) of Theorem 8. In particular, the perturbed Weyl scalars are 
finite at the Cauchy horizon ch + . ■ 

5. The I = 1 Perturbation. 

We return now to the I = 1 perturbation. The treatment is considerably more 
straightforward, but at the loss of full gauge invariance of some of the results. The 
crucial difference for / = 1 is that the metric perturbation quantity kA is no longer 
gauge invariant. We find that under the infinitesmal transoformation generated by 
£ = £,S a dx a , kA transforms as 

k A ^k A -r 2 (r- 2 0,A- 
Furthermore, the equation (|2~%j) no longer holds. However the quantity II is gauge 
invariant, and the equation (|H2*|) holds. Since I = 0, this equation is readily solved 
once La is determined. This is done following the same procedure as for I >2: we find 
La = Ho(v)pkA- It is useful to take a different approach (see the corresponding treatment 
of the odd-parity I = 1 perturbation in PI]). The divergence form of the conservation 
law (j2Uj) indicates the existence of a potential for La- we may write r 2 L^ = €a b J,b- 
Comparison with the previous version of La shows that 7 = j(v) with 7' = — 2X/j,q. The 
advantage of this is that we can now write (J32|) in the form 

e A B (167T 7 -r 4 n) iB = 0, 

yielding r 4 II = 16777(1;) + c, where c is a constant of integration. This form applies 
throughout the spacetime, including the region v < 0, where the background is flat. It is 
appropriate to assume a vanishing matter perturbation ( 7 = 0) in this region, and hence 
the appropriate boundary condition for IT is to take c = 0. Thus the gauge invariant 
perturbation represented by II is completely determined by the matter perturbation 
quantity 7 (or equivalently, /i ). Thus we have 

r 4 II = 16717(1; ), 

and there is no divergence at the Cauchy horizon (except possibly at r = 0, depending 
on the details of 7). 

Consdering the perturbed Weyl scalars, we note that 8^0 and 6^4 vanish identically 
for Z = 1 (as expected: this corresponds to the absence of dipole gravitational radiation). 
5^2 is essentially IT, and so the comments above regarding finiteness apply also to this 
Weyl scalar. From (fKUJ) and (|51jl. it is clear that 5^> \ and 6^3 are not gauge invariant 
for 1 = 1. The effect of the gauge transformation generated by £ is 

^i^^ 1 + 4 ^i(r- 2 0,A/ A , 
r 

^ 3 -^3 + 4^(r- 2 0,An A 
r 

So while we cannot ascribe any direct physical significance to these terms, it is also true 
that any divergence of these quantities is gauge-dependent, and can be removed by a 
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gauge transformation. Thus the situation for I = 1 is the same as that for I > 2: an 
initially finite perturbation remains finite at the Cauchy horizon. 

6. Conclusions 

In this paper, we have studied odd-parity perturbations of self-similar Vaidya spacetime. 
More accurately, the study is of the multipoles of the perturbation, i.e. the coefficients of 
the (scalar, vector, tensor) spherical harmonics with respect to which the perturbation 
quantities may be decomposed. The results are very straightforward to state in rough 
terms: if the perturbation is initially finite, then it remains finite as it impinges on the 
Cauchy horizon. This statement of our results hides the details that are represented by 
Theorems 1-9: the word 'perturbation' refers to the gauge invariant metric, matter 
and Weyl tensor quantities, and 'finite' refers both to integral energy measures and 
pointwise values. Another detail is the meaning of the term 'initial': we slice the 
relevant region of spacetime using hypersurfaces generated by the homothetic Killing 
vector field. These have the advantage of being naturally aligned with the evolution of 
fields on the spacetime in the sense that when we use this slicing and the associated 
time coordinate t, the evolution equations are independent of the space coordinate x. 
The disadvantage is that these slices meet the (singular) scaling origin of the spacetime 
rather than the regular axis. Consequently, one is driven to consider data that vanish 
in a neighbourhood of this point. This is undesirable, as this corresponds to data that 
are supported outside the past light cone: this is clearly not the most general kind of 
data one would like to consider. However, we can circumvent this problem by studying 
a rescaled version of the fields, for example - and most importantly - <fi = r K U. Then 
taking data with <ft e H 3 ' 2 (and one derivative less for its time derivative, along with 
an appropriate specification of the initial matter perturbation), one obtains results for 
which the physical field II does not necessarily vanish at the origin, and for which all 
the finiteness results follow through. 

In previous work, we considered the even parity perturbations of self-similar Vaidya 
spacetime |T2|. Here, it was necessary to take a Mellin transform of the (much more 
complicated) system of perturbation equations. The results for the individual modes 
of the perturbation were as for the perturbations studied here: an initially finite 
perturbation remains finite at the Cauchy horizon. These results and those of the 
present paper provide evidence for the stability of the naked singularity in self-similar 
Vaidya spacetime. As noted in the introduction, this should not however be considered 
a strong challenge to the cosmic censorship hypothesis. Nonetheless, the results do 
indicate the propensity of self-similar naked singularities to survive intact under linear 
preturbations. Our hope is that the approach here can be applied to cases of greater 
physical interest (especially those of perfect fluid [TH] and sigma model [THj collapse) to 
yield insights into cosmic censorship in these cases. 
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